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ABSTRACT 

This thesis develops a block diagram representation of a salient 
pole synchronous generator. This representation in its final form 
provides an accurate means for analyzing the machine by analog com- 
puter methods. 

Linear volt-ampere equations describing the synchronous 
alternator are developed from the basic equations for the various 
machine circuits by means of matrix algebra. Using these linear 
volt-ampere equations and making certain simplifying assumptions, 
a single time constant computer block diagram which represents the 
alternator and its connected load is developed. The nonlinear effects 
of saturation are introduced at a single point in this block diagram. 

The validity of the representation is shown by a comparison of 
steady state and transient phenomena as obtained experimentally from 
the alternator and analytically from the computer representation of 
the alternator. The steady state results show exact-agreement of the 
open circuit condition and only slight divergence in the upper regions 
of saturation for the experimental and computer curves describing 
operation at both zero and unity power factors with rated armature 
current. 

Transient phenomena are measured by application of both load 
steps and field excitation steps for three widely separated conditions 
of loading. Agreement between experimental and computer transients 
is excellent for the load changes and acceptable for the field excita- 
tion changes. 
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Conclusions 



(1) The block diagram developed is a valid representation of the 
synchronous alternator. 

(2) The final representation is in a form readily adaptible to computer 
presentation. 

(3) The experimental and computer steady state characteristics curves 
diverge somewhat at high values of field excitation for both unity 
and zero power factor loads. This divergence will decrease for 
intermediate values of loading. 

Recommendations 

The recommendation is made that the following points should be 

investigated further: 

(1) Methods to reduce the steady state errors at high values of field 
excitation with both zero and unity power factor loads. 

(2) Plausibility of including the effects of armature transients and 
amortisseur subtransients in the alternator block diagram. 

(3) A method to properly account for speed variations in the repre- 
sentation. 

(4) Study of inclusion of unbalanced loads in the generator-load repre- 
sentation. 
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INTRODUCTION 



The transient performance of synchronous machines has been 
the subject of many careful investigations which have resulted in a 
considerable amount of data representing operation under all types of 
conditions. Most of these investigations have been based upon linear 
theory in order to keep the analysis as clear and uncomplicated as 
possible. When included in these analyses, nonlinear effects were 
introduced by adjustment of machine parameters at a particular con- 
dition of operation and were valid only at that condition. 

A more recent investigation (11)(12) by the Energy Conversion 
Group, Massachusetts Institute of Technology, adopted a somewhat 
different approach. This investigation was primarily concerned with 
the voltage regulatory problem of an alternating current generator 
system. Since the various dynamic problems associated with regula- 
tion, such as load changes, faults, and switching, as well as the many 
components of the system itself, are nonlinear, a completely linear 
analysis was impossible. Consequently, one of the purposes of this 
investigation was the development of a block diagram representation 
of the voltage regulatory system, including nonlinear effects, which 
would be suitable for a study of the problem using an analogue computer. 

In the portion of the study concerned with the synchronous 
alternator the Energy Conversion Group was successful in obtaining 
a computer representation of the alternator which proved to be valid 
with a zero power factor load. 



1 



The purpose of this report is to present the steps which were 
followed in the development of a dynamic computer block diagram 
representation of a synchronous alternator, including nonlinear char 
acteristics, which is valid for any type of transient disturbance and 
for loads of power factors from zero to one inclusive. 
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PROCEDURE 



The procedure followed in developing the computer block diagram 
representation of a synchronous alternator and testing its validity can 
be divided into the following steps: 

Development of the Linear Volt- Ampere Equations Describing the 
Synchronous Alternator 

In this portion of the report the electrical equations of motion for 
the stator and rotor were written and, by a series of transformations, 
were changed to direct and quadrature two- axis components (10). The 
equations for the synchronous machine were first developed in motor 
form and were then changed to the generator form. Throughout this 
development matrix algebra was used to provide a compactness of form 
with the transformations following a path from three-phase components 
to Lyon' s symmetrical components, to forward and backward rotating 
components, and finally to Park's direct and quadrature two-axis com- 
ponents. In order to obtain a final set of volt- ampere relationships 
which are linear differential equations with constant coefficients, it 
was necessary to assume constant speed operation. 

Appendix A is a more thorough description of this portion of the 
report. 

Development of a Computer Block Diagram Representation of the 
Synchronous Alternator and its Connected Load 

A block diagram was developed from the volt-ampere equations 
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describing the synchronous alternator. In order to include a load 
representation in this block diagram a vector diagram describing the 
machine in two-axis components was constructed. Since the block 
diagram developed was not suitable for direct application to an analog 
computer, it was necessary to change this block diagram to computer 
form and include a valid representation of the nonlinear characteristics 
of the synchronous alternator. 

Appendix B describes the development of a computer represen- 
tation more completely. 

Experimental Determination of Machine Parameters and Transients 
Phenomena 

After choosing a suitable synchronous generator, the machine 
parameters were determined from the following tests and measure- 
ments: slip test, open circuit test, short circuit test, and a measure- 
ment of the field time constant. This experimental work provided the 
necessary machine parameters which were then converted to per unit 
values in the reciprocal system described by reference 8. Resulting 
from this conversion were the per unit values necessary to complete 
the computer block diagram representation. These are listed below, 

Rj = per unit field resistance including field rheostat 
and external circuit resistances, 

Xmd “ magnetizing reactance, 

= per unit field leakage reactance, 

X^j^ = per unit armature leakage reactance, 

and the nonlinear relation between resultant open circuit field 
current, i^ , and "air gap” field current, i^. 

^oc ag 
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The transient characteristics of the synchronous generator were 
determined at three widely separated load conditions: pure inductive 
load (pF = 0), pure resistive load (pF = 1), and a combination inductive- 
resistive load (pF = 0. 752). The choice of load also determined the 
remaining two blocks, \| + 1 andKj(0), in the load portion of 

the computer representation. For each of these load conditions the 
dynamic effects upon armature terminal voltage, armature current, 
and field voltage were determined by oscillograms when making the 
following setp changes: no load to full load, full load to no load, and 
a field voltage step. 

The final computer block diagram representation was based 
upon the assumption that the machine speed remained constant. During 
the determination of the transient phenomena, the machine speed 
varied upon application of the various step changes, the amount of 
variation depending upon the load condition. An oscillogram was taken 
of this speed transient in order to determine the effect of the speed 
variation upon the other transients measured. 

Appendix C is a more complete description of the experimental 
determination of machine parameters and transient phenomena. 

Computer Determination of Transient Phenomena 

The alternator representation obtained should be suitable for 
analytic study and solution by analog computer methods, and was 
developed with this object in mind. In order to test the validity of the 
alternator representation, the solution of the block diagram should be 
checked with transient and steady- state values obtained by direct ex- 
perimentation on the test machine. An initial computer test was. 
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therefore, to obtain characteristic curves of steady-state relationships 
between field excitation and generator voltage under various conditions 
of loading. These curves should compare with the saturation curves 
experimentally obtained. The computer study of transient phenomena 
was made duplicating the field excitation and load conditions which 
were used in the experimental transient study. 

Appendix D describes the computer analysis more completely. 
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RESULTS 



Figure I presents the block diagram representation of the 
salient pole synchronous generator. This representation was used 
to determine various steady state and transient terminal voltage 
responses from an analog computer study. An experimental study 
was also made duplicating the conditions used for obtaining the com- 
puter responses. This section of the thesis correlates the results 
of the experimental and computer work and proves the validity of 
the developed block diagram. 

Figure XI presents a comparison of the saturation curves 
obtained by both experimental and computer study at no-load and at 
zero and unity-power factors with rated armature current. 

Figures II through X present a comparison of the transient 
responses of the terminal voltage for changes in load and field exci- 
tation as obtained by both computer and experimental studies. The 
initial and final voltages in the case of each transient were adjusted 
on the experimental curves to coincide with the initial and final 
voltages obtained from the computer curves. This was made neces- 
sary because the computer study was carried out assuming constant 
speed, whereas actually the experimental transients were accom- 
panied by generator speed changes. 

The computer representation was developed by neglecting 
armature transients and subtransient effects due to the amortisseur 
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windings. The computer transients reflect this onlission in the 
theoretical development by a sudden change at the point where the 
transient was initiated. These transient and subtransient reactance 
changes are present in the experimental curves near the point of 
transient initiation; however, they die out very rapidly. 

Table I is a comparative tabulation of the time constants 
obtained from the experimental and computer transients as measured 
when neglecting those transients due to armature reactance and 
amortisseur subtransient reactance. 
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TABLE I 



TIME CONSTANT COMPARISON 



Power factor = 0 


Experimental 


Computer 


Full load to no load 


0. 10 


0. 105 


No load to full load 


0. 125 


0. 145 


Field step 


0. 155 


0. 20 


Power factor =1,0 






Full load to no load 


0. 125 


0. 13 


No load to full load 


0. 13 


0. 14 


Field step 


0. 148 


0. 175 


Power factor = 0. 752 






Full load to no load 


0. 11 


0. 10 


No load to full load 


0. 14 


0. 145 


Field step 


0. 145 


0. 195 
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BLOCK DIAGRAM REPRESENTATION FOR THE SALIENT POLE SYNCHRONOUS ALTERNATOR 
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IV 



DISCUSSION OF RESULTS 



Steady State Relationships 

The no load saturation curves obtained by computer and 
experimental methods compared exactly over the entire range of field 
excitation used. This is as might be expected since the circuit at no 
load in the block diagram representation is completely defined by fixed 
quantities, i. e., the nonlinear characteristic, the field voltage, and 
field Impedance. Thus, the no load correlation merely verifies the 
construction of the nonlinear characteristic curve. 

The zero power factor saturation curves for both the computer 
and experimental studies were found to compare exactly as field cur- 
rent was increased until a point was reached slightly above rated ar- 
mature voltage. At that point the computer curve began falling below 
the experimental curve with slowly increasing divergence. The analy- 
tical (computer) representation was developed assuming that the arma- 
ture leakage reactance was a constant quantity equal to the Potier 
reactance (see Appendix C). Thus, the same Potier triangle moved 
parallel to itself would fit between the open circuit and zero power 
factor curves throughout the range of field excitation. 

The zero power factor saturation curve does not have exactly 
the same shape as the open circuit characteristic shifted downward 
and to the right by constant amounts. The principal reason for this 
discrepancy is that the field leakage flux is not a constant under varying 
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conditions of loading. At zero power factor and with higher field 
excitations than those required for rated armature voltage, the field 
current is greater than for the same armature voltage on open cir- 
cuit. This causes a higher field leakage flux and increased saturation 
due to the main air gap flux and the field leakage flux. This effectively 
causes the armature leakage reactance drop to be somewhat less than 
the Potier reactance. In general, the presence of field leakage in a 
salient pole machine causes the open circuit and load characteristic 
curves to draw near each other vertically for high values of saturation. 
This subject is covered with considerable clarity by Fitzgerald and 
Kingsley (see reference 4). 

The correlation between the unity power factor saturation curves 
obtained experimentally and by computer is fairly close, agreeing 
exactly in the lower linear regions and diverging somewhat with in- 
creasing values of saturation. The computer curve gives higher arma- 
ture voltage than the experimental curve at the same field excitation 
when in the saturated region. This discrepancy is largely due to the 
assumption of no saturation of magnetic material in the quadrature 
axis. This assumption is not entirely correct, but because the inter- 
polar air gap spaces have such a high proportion of air in the magnetic 
circuit as compared to that in the direct axis magnetic circuit, the 
value of quadrature axis synchronous reactance, X , was assumed 

M. 

constant and unsaturated. Actually, some degree of saturation does 
exist at higher values of field excitation due to the presence of iron in 
the quadrature axis magnetic circuit, although not as much as in the 
direct axis. The effect of quadi’ature axis saturation is to reduce the 
value of X from its unsaturated value. Since the direct axis terminal 
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voltage, e is directly proportional to X , saturation has the effect 
of reducing e^ and, consequently, the terminal voltage. Thus, the 
experimental saturation curve for unity power factor would rightly be 
expected to be lower than the corresponding computer curve. 

It should be noted that while the computer characteristic for 
unity power factor is higher than the unity power factor experimental 
curve, the reverse is true for zero power factor loads as earlier ex- 
plained. This is evident from a study of figure XI of the results. It 
therefore appears quite likely, that at some intermediate power factor 
between zero and unity, the variation of the analytic steady state values 
from those experimentally determined would be zero. For such a 
power factor the effect of decreasing field reactance with saturation 
would be cancelled by the decreasing value of X^ with saturation. 

It would be possible to reduce the discrepancy between computer 
and experimental curves at unity power factor by inserting a satura- 
tion factor at X^ in the block diagram of figure I of the section on 
results to account for the discrepancy which presently exists. Further 
compensation to account for the small discrepancy between the com- 
puter and experimental zero power factor curves could also be accom- 
plished by insertion of a nonlinear factor at the location of X^j. The 
resulting block diagram would give more accurate analytic results but 
would involve three nonlinear elements which are characteristic of a 
single machine. The relative simplicity of the representation developed 
in this thesis would have to be sacrificed to some degree in order to 
obtain a closer correlation between the experimental and analytic 
results. 
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Transient Relationships 



Figures II through X of the section on results graphically compares 
the transient phenomena observed for the experimental and computer 
studies. The comparisons for full load to no load and for no load to 
full load are quite excellent. The transients in the terminal voltage 
due to a field step are also good, although not agreeing as closely as 
for the load changes. 

The experimental transient is rehlly a three time constant 
curve, one due to field transients, another due to armature transients, 
and the third due to subtransient effects of the amortisseur. The 
subtransient component of the experimental wave is of very short 
duration as compared with that portion caused by the field transient. 
Likewise, the armature transient is of short duration as compared to 
the field transient. The alternator block diagram as developed in this 
thesis has neglected the effects of the subtransient and armature tran- 
sient components (see Appendix B). A study of the transients due to 
load changes in figures II through X shows the discrepancy between 
computer and experimental curves caused by neglecting these com- 
ponents. Their presence can only be noted for a short interval follow- 
ing the initiation of the transients, and they have completely disappeared 
after 0. 1 seconds. After their disappearance, the experimental tran- 
sient closely follows the computer transient. This result justifies 
the assumption that the alternator behaves like a single time constant 
machine. 

In each of the nine transient comparisons the time constant 
associated with the experimental transient was somewhat less than 
the computer time constant. The best explanation for this consistent 
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error is that the measured value of field leakage reactance, Xj.^ was 
probably in error since this was the least reliable parameter measure- 
ment. It was obtained by first measuring the time constant of the field 
transient associated with a step in field voltage with open circuit ar- 
mature terminals. The time constant was then multiplied by the field 
resistance to obtain the field inductance. The leakage term was ob- 
tained by subtracting the mutual inductance term, common to both 
field and armature, from the value of field inductance. Thus, the 
field leakage inductance was the difference between two values close 
in magnitude and, hence, in subject to considerable error. A smaller 
value of field leakage inductance would decrease the computer time 
constant and make its transient approach the experimental transient. 

The discrepancies occurring between computer and experimental 
transients for field steps were greater than for load changes. This 
was probably due to the method used for applying the field steps. In 
order to avoid changing the field resistance, the step was applied by 
shifting the excitation voltage source of 112 volts D. C. to another 
voltage source of approximately 70 volts D.C. A double throw switch 
was employed for this purpose. The switching delay was about 0.02 
seconds during which time the excitation voltage was zero. Thus, 
during the delay period, the armature voltage was dropping toward 
zero as a final asymptote. The resulting experimental transient thus 
dropped more rapidly than would be the case if no switching delay 
were introduced. 

The block diagram developed in this thesis does not account 
for changes in speed. As was previously explained, the initial or 
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final values of the experimental transients were adjusted in each case 
where they were graphically compared with the computer transients. 
In general, any change in load or excitation will cause a generator 
speed variation. These speed changes should be incorporated into 
the block diagram representation in order to extend its employment 
in the study of the voltage regulation problem. A method to properly 
account for speed in the alternator representation should be a subject 
for further study. 

The final representation was developed assuming balanced load 
conditions. Any study of the voltage regulation problem should con- 
sider unbalanced load conditions due to faults in a single phase of the 
generator load. Thus, an additional recommendation for further in- 
vestigation would be to study the possibility of including unbalanced 
loads in the generator-load representation. 
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V 



CONCLUSIONS 



(1) The block diagram development obtained in this thesis is a valid 
representation of the synchronous alternator. This conclusion is 
based on the close correlation obtained betw^een the computer and 
experimental results in both the transient and steady states. 

(Z) The block diagram representation is a comparatively simplified 
form. This simplification results in some slight departures 
from actual steady state values in the saturated regions of the 
magnetization curve. 

(3) At unity pow^er factor, the computer steady state characteristic 
lies above the experimental curve. At zero power factor, the 
computer characteristic lies below the experimental curve. At 
load power factors intermediate between zero and one, the dis- 
crepancy between actual and analytic results would probably be 
less than the maximum error obtained at either 0 P. F. or 1. 0 P. F. 
Since the normal load on a generator is rarely at either of these 
extremes, the resulting correlation of computer and experimental 
steady state data utilizing the alternator representation herein 
presented should be better than that obtained for either unity or 
zero power factors. 
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VI 



RECOMMENDATIONS 



(1) Investigation of methods to reduce the steady state errors at high 
values of field excitation with both zero and unity power factor 
loads. 

(2) Investigation of the plausibility of including the effects of armature 
transients and amortisseur subtransients in the alternator block 
diagram. 

(3) Study of a method to properly account for speed variations in the 
representation. 

(4) Study of inclusion of unbalanced loads in the generator-load repre- 
sentation. 
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APPENDIX A 



DEVELOPMENT OF THE LINEAR VOLT-AMPERE EQUATIONS 
DESCRIBING THE SYNCHRONOUS ALTERNATOR 



Introduction 

In the development of the volt -ampere equations for the 
synchronous alternator the methods outlined in reference 10 have 
been followed. The equations are first developed assuming motor 
action and utilizing the machine structure of figure 1. 

FIGURE (I) 

Physical Structure of Machine 



= (D + 



Figure (I) places the salient field pole structure on the stator 
although in practice the d. c. field is usually placed on the rotor. Since 
relative motion is the only consideration, this representation is valid. 
However, in order to conform to standard synchronous generator prac- 
tice, the equations will later be changed to show generator action and 
the salient field structure on the rotor. 




24 



The development below also makes the following assumptions; 

(1) Geometrical symmetry exists on at least one side of 
the air gap. 

(2) No slot openings exist. 

(3) No saturation occurs. 

(4) Sinusoidal mutual air gap flux density exists. 

(5) Amortisseur windings are neglected. 

Assumption (1) is valid since the machine tested has a cylindrical stator 
structure. Assumptions (2) and (3) arc approximations which will be 
used in this portion of the report in order to develop linear differential 
equations with constant coefficients. Appendix B describes the method 
for the introduction of the nonlinear effects of saturation. Assumption 
(4) is a series approximation which is used even though the rotor of 
the machine tested is a salient member; this approximation can only 
be justified by the results obtained. For the machine tested, assump- 
tion (5) is not valid and is only used in order to simplify the analysis. 

In order to include these windings in the analysis it is necessary to add 
a second winding on the stator member of figure (I). 

The development to be presented also assumes that the machine 
equations can be derived from linear lumped circuits in relative motion. 
Therefore, the three phase belts of the armature are replaced by three 
coils each of N equivalent turns which are sinusoidally distributed in 
space. Matrix algebra is used to aid in the development of the equa- 
tions and reference 10 has been used throughout as a guide. 

Phase Components 

The electrical equations of motion for the salient pole machine 
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of figure (I) in terms of the phase variables are 
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= Po + P2 cos 2(0^ + 0) (9) 
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Finally, in terms of the phase components the equations become 
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Symmetrical Components 



In order to simplfy the equations, the Lyon transformation is used 
to obtain symmetrical components. To accomplish this simplification, 
use is made of the transformations 



a = i 
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a a 
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a a 



where a = e 
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and 
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and the voltages, currents, and impedance are converted to symmetrical 
components as follows: 
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cos(20 + 2d) 
cos 20 
cos(20 + d) 

( 22 ) 

(23) 

(24) 

1 

(25) 

( 26 ) 
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Combining equations (16) through (26) in order to obtain a complete 
volt-ampere equation in Lyon instantaneous symmetrical components, 
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Rotating Components 

In order to obtain linear differential equations with constant 
coefficients the positional dependency must be removed and a speed 
dependency substituted. Then, if constant speed is assumed, the 
equations become linear. For a salient pole machine the necessary 
transformations can only be made by referring the symmetrical mem- 
ber to the salient member: in this instance this will mean that the rotor 
quantities will be referred to the stator. In equation (27) only the +, 

-, and field components are positionally dependent and torque -producing; 
therefore, only these components need to be transformed. In order to 
convert symmetrical components to rotating components the transfor- 
mations shown in equation (28) will be used. 
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With these transformations the voltages, currents, and impendances 
in rotating components become 
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After completing the necessary matrix multiplications the completed 
equations in rotating components become 
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Park's Direct and Quadrature Two -Axis Components 

Since the rotating components are speed dependent with rotor 
quantities referred to the stator, the real direct and quadrature two- 
axis components can be obtained directly from the complex rotating 
components by using the change of variable of equation (32). 
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The volt-ampere equation in two axis components becomes 
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Combining the results of equations (34), (35), and (36), the completed 
equations in Park's two -axis components are 
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Defining the direct axis and quadrature axis inductances as 
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The equations for motor action and the physical structure of 
figure (I) become 
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Generator Equations 

The volt-ampere equations have been derived for motor action 
assuming clockwise rotation of the armature as indicated in figure (I). 
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In order to obtain generator equations the following changes must be 
made: 

(1) Since the field (rotor) rotates in a clockwise 
direction, the armature must revolve counterclock- 
wise or opposite to the direction shown in figure (I). 
Therefore, all speed terms in the matrix must have 
their signs changed. 

(2) The matrix as given is developed in terms of voltage 
drops. In order for the equations to be in terms of 
generator voltage rises, all V terms must have their 
signs changed, i. e. , e = -V. 

(3) The electromagnetic torque for a motor is positive 
while for a generator it is negative. Since the sign 
of the torque is governed by the polarity of the exci- 
tation current, must be made negative for a gen- 
erator. 

Making the changes indicated above, the volt-ampere matrix for a 
generator becomes 
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(40) 



Thus, the linear volt-ampere equations describing the synchronous 
alternator in Park's d-q two-axis components become 

e^ = -(R'^+ pL^)Ij + + Mp Ij (41) 
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= -0 L , I . - (R^ + p L ) I + M 0 1. 
q ad ‘ q q 1 



( 42 ) 



MI^ + (Rj. + pLj)Ij 



(43) 
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APPENDIX B 



DEVELOPMENT OF A COMPUTER BLOCK DIAGRAM 
REPRESENTATION OF THE SYNCHRONOUS ALTERNATOR AND 
ITS CONNECTED LOAD. INCLUDING EFFECTS OF NONLINEARITIES 

Development of Linear Block Diagram 

The linear differential equations derived in Appendix A for a 
salient pole synchronous generator with resistive -inductive loading 
are 




( 1 ) 



e = -0L ,I . - (R + pL )I +M0L 
q dd' q q f 

e^ = -p 2 MI^ + (Rj. + pLj)Ij 
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FIGURE (I) 
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These equations can be represented by a block diagram in which e^ is 
the exciting signal and the direct and quadrature currents and voltages 
feed a load as yet undefined. 

The block diagram obtained above is derived from equations 
which do not account for the subtransient effects due to the amortis- 
seur. These subtransients are of very short duration compared to 
others which are present, and it is often justified to neglect them in 
considering the problems of voltage regulation and system design. 

It should be noted that eq, i^, and i^ are actually direct 
current modulating envelopes of an A. C. carrier. Variation in these 
quantities occur only transientwise. 

The equations representing the load to be applied to the generator 
terminals can be obtained in a straightforward manner from a study of 
the vector diagram describing the volt-ampere relationship of the load 
in the steady state. The steady state relationship is not an exact rela- 
tionship, since any changes in the load signal will cause transients to 
occur. However, these transients, which are unidirectional, do not 
produce appreciable effects and hence are neglected in the same manner 
that the armature transients are neglected as will be shown later. 




I 



'f A 



^ 0 = phase angle between i^ 

and e^ 



6 = power angle between e^ 
and e^ 



Xj^, and are the 

absolute values of load impedance 



FIGURE (II) 
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i , = i sin (6 + 0) 
da 

also: 

e , = e. sin6 = i Z, sin 6 

U t ' . 3. 

• • _ d 

* * ^a ” Zj^ sin 6 



(4d) 



(5d) 

(6d) 



iq = cos (6 + Q) 



= e. cos 6 = i Z, cos 6 
q t a L 



1 = 
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(4q) 

(5q) 

(6q) 



Substituting equation (6) into equation (4), 
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i . = — cos 0 + sin 0 (9d) 

d 2;^ Zj^ 
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'q ~ z. 



i_ = 7^ cos 0 - 7T— sin 0 



(7q) 

(8q) 



(9q) 



Applying equations (9d) and (9q), the load can now be represented 
in block diagram form and added to the block diagram earlier obtained 
for the generator. 
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FIGURE (III) 



Basic Assumptions 

The block diagram represented above does not account for the 
saturation effects normally encountered in a magnetic circuit. To 
account for the saturation of the machine, the nonlinearity is approxi- 
mated by means of a single characteristic which can be placed in such 
a position in the block diagram as to correct for the linear, 
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characteristic block diagram obtained above. 



In addition to neglecting subtransient components, it can also be 
seen that armature transients can be neglected. These transients are 
largely responsible for a second harmonic wave which is not detect- 
able on the experimental transient current wave forms, and also 
account for the direct current transients of the output wave which vary 
from phase to phase depending upon the time in each cycle at which 
the transient signal is introduced. It would thus be quite difficult to 
compare these D. C. transients experimentally and analytically. 



An acceptable procedure is to neglect them since they .do not 

produce appreciable effects upon the average terminal voltage of the 

three phases. This assumption amounts to dropping the pLj, p L , 

Q q 

and the p M terms from the block diagram. 



In a salient pole machine, the reactance of direct axis armature 
reaction, saturates to a much greater extent lhan the corres- 

ponding quadrature axis reactance, X^^^. It is, therefore, a reasonable 
assumption that X^^, and thus L^, does not saturate, but remains 
linear. 



Two final assumptions are that the armature resistance is 
negligible, and that the speed remains constant corresponding to an 
electrical frequency of 60 cps. 



Per Unit System 

In order to simplify the linear block diagram so as to more easily 
handle saturation effects, a per unit system will be used as here 
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defined. The stator quantities based on machine rating are defined as 

BASE ARMATURE QUANTITIES 
e^^ = rated peak phase voltage 
= rated peak phase current 
X = wL = e /i = armature base reactance 

cLO 0.0 3.0 30 




The per unit system as now defined for the rotor circuit is a 
reciprocal system in which the base field current is the field current 
which will induce in each stator phase a voltage equal to i^^. (8) 

The per unit impedances of the elements in the rotor circuit (field 
circuit) are obtained by transferring the corresponding physical im- 
pedances to the stator by an effective stator-rotor turn ratio, and 
dividing the transferred values by This makes X^^ = X^^ = X^^^ 

= X^^ in the per unit system. Thus, the following per unit field 
quantities can be defined: 



i- = that field current which is read directly on the air-gap 
° line of the no-load saturation curve corresponding to 
a per unit armature voltage equal to per unit X^^ 



= 



fo 



^ ®ao ^37T1 ^ 



ao 



Y -IJIL 
^md " 2 L ^ITzl ' 
ao ao 



3 Lf 

x„. = 4 ( 



fo ^2 



mf - 2 '3721^^ 
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3 

R. = 



fo 



f " 2 X. _ '37n ^ 



ao 



ao 



per unit 0 = 1 at rated speed 



p = per unit time differential operator = ^ 



Block Diagram Including Nonlinear Effects 

If the block diagram parameters are expressed in the per unit 
reciprocal system as defined in the preceding section, and the assump- 
tions of Section II are applied, equations (1), (2), and (3) can be written. 



I 

d aq q 



e = X , I- - (X„, + X ,) I , 
q md f al md d 



6f ■ ^^f P^md^^f ■ P^md ^d 



( 10 ) 

( 11 ) 

( 12 ) 



The effect of saturation is introduced in a single characteristic 
which effectively converts an open circuit voltage to a linearized air 
gap voltage. This introduction of the nonlinearity is accomplished by 
changing from the open circuit field current to a hypothetical air gap 
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The complete block diagram which accounts for saturation of the 
magnetizing circuit (assuming X linear) is given below in per unit. 




FIGURE (V) 



The load portion of the block diagram can be further reduced if 
the load parameters are considered constant. For the experimental 
verification of the block diagram, Z, and X are constant, and for a 

J-i cj 

specific load condition, 0 is constant. Thus, e^ = and i^ = 

for constant power factor. 

Thus: 

id = Kj(e)a^ (13) 

and 

e_j=K2(0)e^ (14) 
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where: 



and 



X cos^ 0 

K,{0) = 4 r-r- 

1 Zt (Zt + sin 0 



X'“L 






sin 0 



K^iO) = 



X cos 0 
3 

Z, + sin 0{X ) 
L q' 



The per phase terminal voltage e^ can be written 










( 15 ) 



( 16 ) 



(17) 



Thus, the generator and load can be represented as 




FIGURE (VI) 



Computer Representation 

The block diagram of figure (VI) is not directly solvable by 
analogue computer techniques because of the presence of differentiators. 
The representation of a differentiation on an analogue computer creates 
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an unusual amount of "noise” in the solution. It is therefore customary 
to transform a diagram involving differentiators into an equivalent 
diagram employing integrators only. The computer representation 
obtained is given in figure (VII). It should be noted that the saturation 
characteristic used is the inverse of that depicted in figure (VI). 




FIGURE (VII) 
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APPENDIX C 



EXPERIMENTAL DETERMINATION OF MACHINE PARAMETERS 
AND TRANSIENT PHENOMENA 



Machine Parameters 

The various machine parameters were determined from the 
following tests and measurements: slip test, open circuit test, short 
circuit test, and a measurement of the field time constant. An oscillo- 
graph was used to record both the slip test and field time constant 
measurements. As shown by the calculations below, the direct axis 
synchronous reactance, X^, was found from the short circuit and open 
circuit tests (see figure II). The quadrature axis synchronous reactance, 
Xq, was then calculated from the value of X^ already found and from the 
ratio Xq/X^ as determined from the slip test (see figure I). 




FIGURE I. SLIP TEST 
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LINE TO LINE O.C VOLTAGE C»«*) 
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From the short circuit and open circuit tests. 



= = 246 

43,2/7 



3. Z85 ohms. 



( 1 ) 



From the slip test. 



\ _ 1.04/1.6 

Xd 1.06/1.08 



0. 662 



Therefore 



(2) 



Xq = 0. 662 = (0. 662){3. 285) = 2. 175 ohms. (3) 

The direct and quadrature synchronous reactances, X^ and X^, 
can be expressed in per unit as the sum of an armature leakage reac- 
tance X^j^ and a reactance X^^ corresponding to the direct axis arma- 
ture reaction flux: 



^d = X^l ^ad 



X^ = x^, + X^^ 
q ai aq 



(4) 

(5) 



In a similar manner the per unit field self reactance X|.^ can also be 

expressed as the sum of a field leakage reactance X^^ and the field 

magnetizing reactance which in the reciprocal per unit system 

(8) is equal to X ,. 

' md 



Xff = Xj^ + Xjj^ - X^^ + X^d 



( 6 ) 



In order to obtain an adequate block diagram representation it 

was necessary to determine the alternator field resistance R^, a proper 

load representation, a method for insertions into the block diagram the 

nonlinear effects of saturation, and values for X,,, X ,, and X„ , which 

li al md 
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have already been defined. For the computation of the values above, 
a reciprocal per unit system (8) was used in which base field current 



is so defined that the reactance of direct axis armature reaction, 

is the equivalent per unit magnetizing reactance between armature and 
field, Xj^j. Therefore, base field current can be read directly on the 
air gap line of the no-load saturation curve corresponding to a per unit 
armature voltage equal to per unit X^^. 

Of the parameters needed for the block diagram, Appendix B 
describes the methods used to determine a load representation and the 
nonlinear saturation characteristic. The calculations below determine 
the remaining parameters: R^, X^j^, and X^^. In all cases per 

unit values are found by the reciprocal system described in reference 
8 and are referred to the armature. 



The magnetizing reactance X^^ and the armature leakage 
reactance X^j are found from the Potier triangle of figure (III). The 
distance ac is proportional to X^, be is proportional to = X^^, and 
ab is proportional to It follows that 



= 0. 7895 

X^j ac 



( 7 ) 



and 

= ^ = 0. 2105 
Xd ac 



(8) 



Using the value for Xd already determined. 



Xad = 0. 7895 X^ = 2. 59 ohms 


(9) 


X^j = 0. 2105 Xd = 0. 692 ohms 


(10) 



49 



LINE TO LINE ARMATURE VOLTAGE - RMS 




and in per unit values 



X 



X ^ = X 



ad 



md ad x 



= 0. 592 per unit 



( 11 ) 



ao 



X 



‘■al 



al X, 



= 0.1 576 per unit 



( 12 ) 



ao 



where the base armature reactance, X = e /i = 4. 39 ohms and 

ao ao ao 

e and i are nameplate values (peak), 
ao ao ^ \t' I 

The field leakage reactance X^j^ was obtained by finding the field 

time constant from the application of a step change in field excitation 

and a record of the resulting transient on an oscillogram (see figure IV). 




FIGURE (IV) 
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Thus, knowing the field time constant and having measured the field 
resistance Rj., the field inductance could be calculated. 

Lf = Tf Rf = 6. 96 henries (13) 

The determination of per unit field self- reactance follows. 

= I = 0. 69 per unit (14) 

^ao ao 



where L = X„ /ca 
ao ao 

i = peak armature rated current 
ao ^ 

i- = base field current read on the 

air gap line of the no-load satura- 
tion curve corresponding to a per 
unit armature voltage equal to per 
unit X^j. 

= 1.7 amps 



Since Xjj Xjj + 



Xfi = 0. 69 - 0. 592 = 0. 098 per unit 



(15) 



The only per unit values remaining to be found for a complete 
block diagram representation can be calculated from the following 
formulas: 



Xd 

Xj = = 0. 749 per unit 

(per unit) ^ao 



X 



X 



^(per unit) 



*^(ohms) 

~1T 



0. 495 per unit 



ao 



( 16 ) 



(17) 
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R 



R 



= 2. 63 X lO"'^ Rj 



(ohms) 



(18) 




(19) 



Transient Phenomena 

The dynamic performance of the synchronous alternator was 
tested at three widely varying load conditions: pure inductive load 
(pF = 0), pure resistive load (pF = 1)> and a combination inductive - 
resistive load (pF = 0. 752). At each load condition step changes were 
applied to the alternator and the transient effects upon armature ter- 
minal voltage, armature current, and field voltage were recorded by 
means of oscillograms. The step changes consisted of going from full 
load to no load, going from no load to full load, and applying a step in 
the field excitation. 

The oscillograms of figure (V) through (XIII) indicate the steady 
state values before and after each step change and the time constants 
associated with the transient phenomena for the three load conditions 
used. 
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Experimental Transients - pF = 0 




Figure V 

Full Load to No Load 
e. = 0. 988 per unit 

h 

e. = 1 . 372 per unit 
^2 

T = 0. 10 sec 




Figure VI 

No Load to Full Load 

e. = 1 . 363 per unit 
1 

e, = 0. 982 per unit 
^2 

T = 0. 125 sec. 




Figure VII 

Field Step 

e, =1.0 per unit 

e. = 0. 705 per unit 
^2 

T = 0. 1 55 sec. 
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Experimental Transients - PF =1.0 








Figure VIII 

Full Load to No Load 

e. =1.0 per unit 
1 

e. = 1 . 327 per unit 
^2 

T = 0. 125 sec. 




Figure IX 

No Load to Full Load 

e, = 1 . 325 per unit 
1 

e. = 0. 988 per unit 
^2 

T = 0. 130 sec. 
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Figure X 
Field Step 
e. =1.0 per unit 
e. = 0. 734 per unit 

h 

T = 0. 148 sec. 



Experimental Transients - PF = 0. 752 
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Figure XI 

Full Load to No Load 
e = 0. 991 per unit 

h 

e = 1 . 446 per unit 
^2 

T = 0.11 sec. 
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Figure XII 

No Load to Full Load 

e. = 1 . 423 per unit 
1 

e. = 0. 978 per unit 
2 

T = 0. 14 sec. 




Figure XIII 

Field Step 

e, = 1 . 009 per unit 

e. = 0. 737 per unit 
^2 

T = 0. 145 sec. 
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APPENDIX D 



COMPUTER DETERMINATION OF TRANSIENT PHENOMENA 



In order to check the validity of the synchronous generator 
representation developed in Appendix B, it is necessary to carry out 
both experimental and analytic tests which can be correlated in the 
transient and steady states. The purpose of this appendix is to obtain 
the analytic results of various changes in excitation and load conditions 
on the generator output. 

One of the most difficult problems in the analysis and synthesis 
of the synchronous generator arises from consideration of saturation 
and other nonlinearities. The method employed in this thesis consists 
of concentrating all the generator nonlinearities at a single point in a 
block diagram which can be represented on an analog computer. 

The principle of operation of the analog computer is based on the 
fact that analogies can be drawn among the various types of physical 
systems when the differential equations describing the systems are 
comparable. Of the various types of analog computers, the one most 
frequently used is the d-c analog computer in which direct voltages are 
used to represent the dependent variables. The simple computing net- 
works in the d-c analog computer are made up of RLC networks and 
d-c amplifiers. 

Although the analog computer can carry out the differentiation 
operation, the use of differentiators is avoided when possible because 
they tend to amplify signal noise and cause unwanted high frequency 
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oscillations. Usually when differentiators are present in a circuit 
they can be replaced with integrators by proper adjustment of the 
analog circuit. 

There are several ways of representing system nonlinearities on 
the analog computer. The one resorted to in the solution of the syn- 
chronous generator problem is to construct a nonlinear characteristic 
curve relating the field current (1^^,) which causes an open circuit 
armature voltage, to a field current, which generates the same 

armature voltage on the air gap line of the open circuit characteristic 
curve. From the curve obtained, a photographic plate is prepared for 
use in the computer function generator. A signal input to the function 
generator causes a cathode ray beam to ride on the photographic con- 
tour generating a continuous output signal. Figure (1) shows the photo- 
graphic plate used in the computer work described in this appendix. 




FIGURE (I). COMPUTER FUNCTION GENERATOR PLATE 
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The block diagram of figure (VII), appendix B, must be 
represented on the analog computer in terms of the computer tools 
available; i.e., integrators, summers, potentiometers, voltages, etc. 
This representation is given in figure (II). 

A shift in the time scale of the integrators was performed by 
dividing all integrator input potentiometer settings by the time scale, 

Q^. Thus, if were 100, the computer time to perform a given opera- 
tion would be one hundred times as long as for the real time alternator 
operation. This shift in time scale was made necessary for two reasons 
to allow the transients to be recorded at a sufficiently slow speed that 
the dynamic response of the recording instruments would not be ex- 
ceeded, and also to adjust the signal level of the integrator output to 
insure that the maximum signal level of 100 volts would not be exceeded. 

The differential operator for this per unit work is taken equal to 
d/d(w^t) where Wq is the electrical frequency of the generator, 120 •77' 
radians per second. Thus, p = (l/l207T)p' where p' is d/dt. The 
integration term, 1 /p, therefore equals (120 TT)(I /p' ). In order to 
reduce the per unit block diagram integration to a pure integration, the 
integrator output signal is then divided by 120'/7T The RC network 
shown at the input No. 8 summer is a stabilization network used to 
prevent system oscillation. 

Transients are initiated by throwing a switch on the computer. 
The switch movements cause a step in the field voltage or a change in 
the load condition between full load and no load. 

A preliminary step in the use of the computer is to clieck the 
steady- state values of field current vs. armature voltage for various 
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FIGURE n. COMPUTER CIRCUIT DIAGRAM 



conditions of loading. This step consisted of obtaining the generator 
saturation curves for no load, 1. 0 power factor, and zero power fac- 
tor. These curves are given in figure (III). 

The conditions under which the experimental work of appendix C 
was determined were duplicated on the computer. The computer poten- 
tiometer settings which were used for each solution are tabulated in 
table (I). The transient solutions obtained, with their steady-state and 
transient values, are pictured in figures (IV) through (XII). 
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Computer Transients - PF = 0 




Figure IV F\ill Load to No Load 

e, = 1.205per unit t = 0. 105 sec. 

e = 1 . 332 per unit 
^2 




Figure V No Load to Full Load 



e. =1. lOB per unit t = 0. 145 sec. 

^1 

e. = 0. 995 per unit 

h 




Figure VI Field Step 

e. = 1. 00 per unit t = 0. 20 sec. 

'l 

e. = 0. 66 per unit 
^2 
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Computer Transients - PF =1.0 




Figure VII Full Load to No Load 



e, =1.0 per unit t = 0.13 sec. 

1 

e. =1.11 per unit 
^2 




Figure VIII No Load to Full Load 

e = 1. 12 per unit t = 0. 14 sec. 

1 

e. =1.01 per unit 
^2 







Figure IX Field Step 

e, = 1 . 00 per unit t = 0. 175 sec. 

1 

e. = 0. 62 per unit 
2 
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Computer Transients - PF = 0. 752 




e. = 1. 142 per unit t = 0. 10 sec 

h 

e. = 1. 28 per unit 
2 




e. =1.0 per unit 
^2 
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Figure XII Field Step 



e. = 1. 00 per unit t = 0. 195 sec. 

1 

e. = 0. 60 per unit 
2 
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STEADY STATE POTENTIOMETER SETTINGS 
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STEADY STATE POTENTIOMETER SETTINGS (CONT* D) 




I 
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TABLE I (CQNT' D) 

TRANSIENT POTENTIOMETER SETTINGS 
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TRANSIENT POTENTIOMETER SETTINGS (CONT* D) 
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APPENDIX E 



SUMMARY OF DATA AND CALCULATIONS 



Definition of Parameters 




X 



ff 



Tm 



X 



fl 



X. 



X, 



ao 



X 



al 



X 



ad 



aq 



X 



md 



M 
L , 



= field voltage 

= direct axis armature voltage 
= quadrature axis armature voltage 
= peak armature rated voltage, line to line 
= || “ armature terminal voltage 

= field current 

= direct axis armature current 

= quadrature axis armature current 

= peak armature rated current, per phase 

= ^ i^ + i^ = armature terminal current, per phase 

= field self reactance 

= field magnetizing reactance 

= field leakage reactance 

= direct axis synchronous reactance 

= quadrature axis synchronous reactance 

= e„^/l = base armature reactance 

ao ao 

= armature leakage reactance 
= reactance of the direct axis armature reaction 
= reactance of the quadrature axis armature reaction 
= magnetizing reactance between armature and field 
= mutual inductance between armature and field 
= direct axis synchronous inductance 
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r 

II 


quadrature axis synchronous 


inductance 


Rf = 


field resistance 








open circuit field time constant 




- 


load impedance 








rated electrical frequency 






Summary 








Parameter Actual values • 




Per unit values 


®ao 


220 volts 




■ 


'ao 


28. 9 IfT amps 






^ao 


4. 39 ohms 






l^ao 


0.01164 henries 






*^0 


120 ir rad/sec 




1.0 




112 volts 




0. 0173 




3. 285 ohms 




0. 749 




2. 175 ohms 




0.495 


^ad 


2. 59 ohms 




0. 592 


^md 






0. 592 


^al 


0 . 692 ohms 




0. 1576 






X 



fl 



0. 69 
0.098 



0.0417 sec 
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APPENDIX F 



EQUIPMENT 



Generator 

Westinghouse synchronous alternator, salient pole (rotor), 
with amortisseur windings. 



KVA 


11 


% pF 


90 


Volts 


220 


Excitation Volts 


125 


Amps 


28. 9 


Excitation Amps 


5 


Cycles 


60 


Style No. 


897363 


RPM 


1200 


Serial No. 


l/5c3982 



Driving Motor 

Westinghouse direct current motor, 

Volts 230 

Open Intermittent H. P. 15 
Speed 550/1650 

Computer 

Reeves Electronic Analog Computer (REAC). 



shunt wound. 

Serial No. 955771 
Type RF 10 

Form A17 
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